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Abstract 

A comprehensive graph theoretical and finite geometrical study of the commutation relations between the 
generalized Pauli operators of A-qudits is performed in which vertices/points correspond to the operators 
and edges/lines join commuting pairs of them. As per two-qubits, all basic properties and partitionings of 
the corresponding Pauli graph are embodied in the geometry of the generalized quadrangle of order two. 
Here, one identifies the operators with the points of the quadrangle and groups of maximally commuting 
subsets of the operators with the lines of the quadrangle. The three basic partitionings are (a) a pencil 
of lines and a cube, (b) a Mermin's array and a bipartite-part and (c) a maximum independent set and 
the Petersen graph. These factorizations stem naturally from the existence of three distinct geometric 
hyperplanes of the quadrangle, namely a set of points collinear with a given point, a grid and an ovoid, 
which answer to three distinguished subsets of the Pauli graph, namely a set of six operators commuting 
with a given one, a Mermin's square, and set of five mutually non-commuting operators, respectively. The 
generalized Pauli graph for multiple qubits is found to follow from symplectic polar spaces of order two, 
where maximal totally isotropic subspaces stand for maximal subsets of mutually commuting operators. 
The substructure of the (strongly regular) A-qubit Pauli graph is shown to be pseudo-geometric, i.e., 
isomorphic to a graph of a partial geometry. Finally, the (not strongly regular) Pauli graph of a two-qutrit 
system is introduced; here it turns out more convenient to deal with its dual in order to see all the parallels 
with the two-qubit case and its surmised relation with the generalized quadrangle Q(4, 3), the dual of VF(3). 

PACS Numbers: 03.67.-a, 03.65.Fd, 02.10.Hh, 02.40.Dr 
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Symplectic Polar Spaces - Finite Projective (Ring) Geometries 



1 Introduction 



The intricate structure of commuting/non-commuting relations between A-qubit observables may serve as 
a nice illustration of the distinction between the quantum and the classical and failure of classical ideas 
about measurements. A deeper understanding of this structure is central to the explanation of quantum 
peculiarities such as quantum complementarity, quantum entanglement as well as other related concep- 
tual (or practical) issues like no-cloning, quantum teleportation, quantum cryptography and quantum 
computing, to mention a few. Many "strange features" of finite quantum mechanics are linked with two 
important open theoretical questions: finding complete sets of mutually unbiased bases [T] and/or solving 
the Kochen-Specker theorem in relevant dimensions [2]. Both problems are tricky and difficult due to a 
large number of the observables involved. Already for a two-qubit system, there are as many as fifteen 
operators — tensor products of the four Pauli matrices. This set can be viewed as a graph if one regards 
the operators as vertices and joins any pair of commuting ones by an edge. The two-qubit Pauli graph, 
henceforth referred to as V[2,2], is regular of degree six, that is, every observable commutes with other 
six; one of its subgraphs, frequently termed as a Mermin's square, has already been thoroughly studied 
due to its relevance to a number of quantum "paradoxes" [2J[3]. For A-qubits (A-qutrits), A > 2, the 
corresponding graphs V[2, N] (P[3, A]) are endowed with 4^ — 1 (9 N — 1) vertices. One of their partitions 
features 2+1 (3 N + 1) maximally commuting sets of 2 N — 1 (3^ — 1) operators each and is intimately 
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related to the derivation of the maximum sets of mutually unbiased bases in the corresponding dimensions 

This paper aims at an in-depth understanding of the properties of the ./V-qudit Pauli graphs by employ- 
ing a number of novel graph theoretical and finite geometrical tools. It is organized as follows. Sec. first 
lists basic notions and definitions of graph theory and then introduces the relevant finite geometries. The 
latter start with the ubiquitous Fano plane, continue with other remarkable finite projective configurations 
(e. g., Pappus and Desargues) and related subspaces, and ends with more abstract and involved structures, 
such as generalized polygons and (symplectic) polar spaces. Sec. [3] introduces the two-qubit Pauli graph 
and discusses its basic properties. The graph's three basic factorizations are then examined in very detail 
and their algebraic geometrical origin is pointed out: first, in terms of the three kinds of the geometric 
hyperplanes of the generalized quadrangle of order two, second in terms of the projective lines over the 
rings of order four and characteristic two residing in the projective line over Z% x2 [6]. Sec. [4] discusses a 
self-similarity of the 7V-qubit graph; one shows that its structure is that of the symplectic polar spaces of 
order two [7] and strongly regular graphs associated with them. Finally Sec. [5] deals with some properties 
of the two-qutrit Pauli graph V[3, 2] and muses about possible finite geometry behind it. 

2 Graphs and geometry 
2.1 Excerpts from graph theory 

A graph G consists of two sets, a non-empty set V{G) of vertices and a set E(G) of two element subsets of 
V(G) called edges, the latter regarded as joins of two vertices. Alternatively, vertices are also called points 
and edges also lines [H [9l [10]. Two distinct vertices of G are called adjacent if there is an edge joining 
them; similarly, two distinct edges with a common vertex are called adjacent. If one vertex belongs to one 
edge both are said to be incident. The adjacency matrix A — [ay] of a graph G with |V(G)| = v vertices 
is an v x v matrix in which ay = 1 if the vertex Vi is adjacent to the vertex Vj and ay = otherwise. The 
degree D of a vertex in a graph G is the number of edges incident with it; a regular graph is a graph where 
each vertex has the same degree. A strongly regular graph is a regular graph in which any two adjacent 
vertices are both adjacent to a constant number of vertices, and any two non adjacent vertices are also 
both adjacent to a constant, though usually different, number of vertices. The graph spectrum spec{G) is 
composed of the eigenvalues (with properly counted multiplicities) of its adjacency matrix. For a regular 
graph, the largest eigenvalue equals the degree of the graph and the absolute value of any other eigenvalue 
is less than D. 

A subgraph of G is a graph having all of its vertices and edges in G. For any set S of vertices of G, the 
induced subgraph, denoted (S), is the maximal subgraph G with the vertex set S. A vertex and an edge 
are said to cover each other if they are incident. A set of vertices which cover all the edges of a graph G 
is called a vertex cover of G, and the one with the smallest cardinality is called a minimum vertex cover. 
The latter induces a natural subgraph G' of G composed of the vertices of the minimum vertex cover and 
the edges joining them in the original graph. An independent set (or coclique) / of a graph G is a subset 
of vertices such that no two vertices represent an edge of G. Given the minimum vertex cover of G and 
the induced subgraph G", a maximum independent set I is defined from all vertices not in G' . The set G' 
together with / partition the graph G. 

Two graphs G and H are isomorphic (written G = H) if there exists a one-to-one correspondence 
between their vertex sets which preserves adjacency. An invariant of a graph G is a number associated 
with G which has the same value for any graph isomorphic to G. A complete set of invariants would 
determine a graph up to isomorphism, yet no such set is known for any graph. The most important 
invariants for a graph G are the number of its vertices v = \V{G)\, the number of its edges e = \E(G)\, the 
degree at each vertex, its girth g(G), i. e., the length of a shortest cycle (if any) in G, its diameter and its 
(vertex) chromatic number. The distance between two points in G is the length of the shortest path joining 
them, if any. In a connected graph, distance is a metric. A shortest path is called a geodesic and the 
diameter of a connected graph is the length of the longest geodesic. A coloring of a graph is an assignment 
of colors to its points so that no two adjacent points have the same color. A c-coloring of a graph G uses 
c colors. The chromatic number k(G) is defined as the minimum c for which G has a c-coloring. 

Quite often the structure of a given graph can be expressed in a compact form, in terms of smaller 
graphs and operations on them. Graph union, graph product, graph composition and graph complement 
are a few [§]. The complement G of a graph G has V(G) as its vertex set, and two vertices are adjacent in 
G if they are not in G. We will also need the concept of the line graph L(G) of a graph G, i. e., the graph 
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which has a vertex associated with each edge of G and an edge if and only if the two edges of G share a 
common vertex. 

2.2 Graphs and finite geometries 

A finite geometry may be defined as a finite space S — {P, L} of points P and lines L such that certain 
conditions, or axioms, are satisfied [TT]. One of the simplest set of axioms are those defining the so-called 
Fano plane: (i) there are seven points and seven lines, (ii) each line has three points and (iii) each point is 
on three lines. The Fano plane is a member of several communities, some of them of great relevance to the 
structure of an iV-qubit system. It is, first of all, a near linear space, that is a space such that any line has 
at least two points and two points are on at most one line. The Fano plane is also a linear space for which 
the second axiom "at most" can be replaced by "exactly" . More generally, a projective plane is a linear 
space in which any two lines meet and there exists a set of four points no three of which lie on a line. The 
projective plane axioms are dual in the sense that they also hold by switching the role of points and lines. 
In a projective plane every point/line is incident with the same number k + 1 of lines/points, where k is 
called the order of the plane. It has been long conjectured that a projective plane exists if and only k is 
a power of a prime number and this conjecture was related to the existence of complete sets of mutually 
unbiased bases for iV-qudits [T^]. The Fano plane is, in fact, the smallest projective plane, having order 
k = 2. Projective planes of order k can be constructed from 3-dimcnsional vector spaces over finite fields 
F^; such planes are necessarily Desarguesian, but there also exists non-Desarguesian planes which do not 
admit such a coordinatization. 

The Fano plane belongs also to a large family of projective configurations, which consist of a finite set of 
points and a finite set of lines such that each point is incident with the same number of lines and each line 
is incident with the same number of points. Such a configuration may be denoted (v a , e^), where v stands 
for the number of points, e for the number of lines, a is the number of lines per point and b the number of 
points per line. If the number of points equals the number of lines one simply denotes a configuration as 
(v a ), although it is not, in general, unique. A configuration is said to be self-dual if its axioms remain the 
same by interchanging the role of points and lines. The Fano plane is a configuration (73). We will soon 
meet other two distinguished projective configurations: the Pappus configuration (93) and the Desargues 
configuration (IO3). All the three configurations are self-dual. Any configuration may also be seen as a 
regular graph by regarding its points as vertices and its lines as edges. 

Recently, another class of finite geometries was found out to be of great relevance for two-qubits - 
projective lines defined over finite rings instead of fields O [13l HU [15] . Given an associative ring R with 
unity and GL(2,R), the general linear group of invertible two-by-two matrices with entries in R, a pair 

(a, (3) is called admissible over R if there exist j,6 £ R such that f a ^ ^ 6 GL2(i?.). The projective 

line over R is defined as the set of equivalence classes of ordered pairs (ga,g(3), where g is a unit of R 
and (a,/3) admissible |16[ I17j. Such a line carries two non-trivial, mutually complementary relations of 
neighbor and distant. In particular, its two distinct points X: (get, g(3) and Y: (gj, gS) are called neighbor 

if { 01 ^ j (£ GL2(i?) and distant otherwise. The corresponding graph takes the points as vertices and 

its edges link any two mutually neighbor points. For R = (the graph of) the projective line lacks any 
edge, being an independent set of cardinality k + 1, or a (k + l)-coclique. Edges appear only for a line 
over a ring featuring zero-divisors, and their number is proportional to the number of zero-divisors and/or 
maximal ideals of the ring concerned (see, e.g., [I3]-[T7] for a comprehensive account of the structure of 
finite projective ring lines). Projective lines of importance for our model will be, as already mentioned in 
Sec. 1, the line defined over the (non-commutative) ring of full 2x2 matrices with coefficients in Z2, as 
well as the lines defined over three distinct types of rings of order four and characteristic two [6]. 

A linear space such that any two-dimensional subspace of it is a projective plane is called a projective 
space. The smallest non trivial exemple (other than the Fano plane) is the binary three dimensional space 
PG(3, 2) of which two-dimensional subspaces are Fano planes. A generalized quadrangle is a near linear 
space such that given a line L and a point P not on the line, there is exactly one line K through P that 
intersects L (in some point Q) [18]. A finite generalized quadrangle is said to be of order (s, t) if every line 
contains s + 1 points and every point is in exactly t+1 lines and it is called thick if both s > 1 and t > 1; 
otherwise, it is called slim. If s = t, we simply speak of a quadrangle of order s. A generalized quadrangle 
of order (s, 1) or (l,t) is called a grid or a dual grid, both being slim. The simplest thick generalized 
quadrangle, usually denoted as W(2), is of order 2; it is a self-dual object featuring 15 points/lines and a 
cornerstone of our model. 
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Table 1: The product properties between any two Pauli operators of two-qubits; = I2. 



Further concepts closely related to a projective space are those of a subspace and of a geometric 
hyperplane. A set of points in a projective space is a subspace if and only if for any line L the set 
contains no point, one point, or all the points of L. More restrictively, a geometric hyperplane H of a 
finite geometry is a set of points such that every line of the geometry either contains exactly one point of 
-ff , or is completely contained in H. 

Last but not least, we need to introduce the concept of a polar space. A polar space S = {P, L} is 
a near-linear space such that for every point P not on a line L, the number of points of L joined to P 
by a line equals either one (as for a generalized quadrangle) or to the total number of points of the line 
A polar space of rank N (N > 2) can also be seen [TO] as a set {P} of points, together with certain 
subsets, called subspaces, such that: (a) every subspace, together with its own subspaces, is isomorphic to 
the projective space PG(d, q) over the finite field F q and of dimension d at most N—l, (b) the intersection 
of two subspaces is a subspace, (c) for each point P not in a subspace R of dimension N—l, there is a 
unique subspace S of dimension N—l such that R n S is (N — 2)-dimensional, and (d) there are at least 
two disjoint subspaces of dimension N—l. A polar space of rank two is a generalized quadrangle. A 
particular class of higher-rank, N > 2, polar spaces called symplectic polar spaces are, as already outlined 
in [7], the geometries behind (strongly regular) multiple-qubit Pauli graphs (Sec. [4]). 



3 The Pauli graph of two-qubits 

Let us consider the fifteen tensor products u% ®cjj, i,j € {1,2,3,4} and ^ (1,1), of Pauli matrices 
<?i = {h,o- x , o-y, a z ), where 7 2 = ( \ J ] , tr a = ( ? n)' cr ^ = (n ) and °v = ia x^z, label them 



as follows 1 = h^cfx, 2 = h®o- y , 3 = h®o- z , a = <J x (g)I 2 , 4 = a x ®a x . ..,b = €r y ®h,. . . , c = a z ®I 2 ,. ■ ■ , 
and find the product and the commutation properties of any two of them — as given in Table 1 and Table 
2, respectively. Joining two distinct mutually commuting operators by an edge, one obtains the Pauli graph 
"P[2,2] with incidence matrix as shown in Table 2. After removing the triple {a,b,c} of the "reference" 
observables, the incidence matrix can be cast into a remarkably compact form (Table 3) which makes use 
of three 3x3 matrices: O (the "zero" matrix), A (the identity matrix) and A (the matrix complementary 
to A). The main invariants of V[2,2] and those of some of its most important subgraphs are listed in 
Table 4. As it readily follows from Tables 1-3, V[2,2] is 6-regular and, so, intricately connected with the 
complete graphs K n , n = 5, 6 or 7. First, one checks that V[2,2] = L(Kq), i.e., it is isomorphic to the 
complement of line graph of Kq. Next, computing its minimum vertex cover (Table 4), one recovers the 
Petersen graph PG = L(K§). Finally, V[2, 2] is also found to be isomorphic to the minimum vertex cover 
of L(Kj). Now, we turn to remarkable partitionings/factorizations and the corresponding distinguished 
subgraphs of V[2,2]. 
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Table 2: The commutation relations between pairs of Pauli operators of two-qubits aka the incidence 
matrix of the Pauli graph V[2,2]. The symbol "0"/"l" stands for non-commuting/commuting; although 
the diagonal should feature l's (every operator commutes with itself), we put there 0's for the reason which 
will become apparent from the text. 
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Table 3: Structure of the incidence matrix of "P[2, 2] after removal of the triple of operators {a, b, c}. 



3.1 The "Fano pencil" FP and the cube CB 

We shall first tackle the 7+8 partitioning of the graph which can, for example, be realized by the following 
subgraphs/subsets: FP = (1, 2, 3, a, 4, 5, 6) and CB = (b, 7, 8, 9, c, 10, 11, 12). The subgraph FP can also 
be regarded as a line pencil in the Fano plane [3l [20] as well as a hyperplane of W(2) [6]; the number 
of choices for this partitioning is obviously equal to the number of the vertices of the full graph (see [3] 
for another choice). A CB is also the generalized Petersen graph G(4, 1) [SJ. Employing Table 1, it is 
easy to observe that two vertices on one line of FP map to the third one on the same line, i.e., l.a = 4, 
2. a = 5 and 3. a = 6. The three observables are found to share a common base of 4-dimcnsional vectors; 
for this particular choice, the lines in the Fano pencil FP feature unentangled 2-qubit bases. In addition, 
an edge of CB is mapped to a vertex of FP, e. g., 8.10 = 6, 8.12 = —4, etc. In particular there is a closed 
path of length 6 (shown with thick lines) in the cube graph CB which features six bases of entangled 
states. It is worth mentioning here that in [3] the projective lines over direct product of rings of the type 
Z 2 XTl , n = 2,3,4, were used to tackle this kind of partitioning. With these lines it was possible to grasp 
the structure of the two subsets, but not the coupling between them; to get a complete picture required 
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Table 4: The main invariants of the Pauli graph V[2, 2] and its subgraphs, including its minimum vertex 
covering MVC isomorphic to the Petersen graph PG. For the remaining symbols, see the text. 
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Figure 1: Partitioning of V[2, 2] into a pencil of lines in the Fano plane (FP) and a cube (CB). In FP 
any two observables on a line map to the third one on the same line. In CB two vertices joined by an 
edge map to points/ vertices in FP. The map is explicitly given for an entangled path by labels on the 
corresponding edges. 
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Figure 2: Partitioning of "P[2,2] into an unentangled bipartite graph (BP) and a fully entangled Mermin 
square (MS). In BP two vertices on any edge map to a point in MS (see the labels of the edges on a 
selected closed path). In MS any two vertices on a line map to the third one. Operators on all six lines 
carry a base of entangled states. The graph is polarized, i.e., the product of three observables in a row is 
— 14, while in a column it is +I4. 



employing a more abstract projective line with a more involved structure [6]. 
3.2 The Mermin square MS and the bipartite part BP 

We shall focus next on the 9+6 partitioning which can be illustrated, for example, by the subgraphs 
BP = (1, 2, 3, a, b, c) and MS = (4, 5, 6, 7, 8, 9, 10, 11, 12). The BP part is easily recognized as the bipartite 
graph K[3, 3], while the MS part is a 4-regular graph. There is a map from the edges of BP to the vertices 
of MS, and a map from two vertices of a line in MS to the third vertex on the same line. The bases 
defined by two commuting operators in BP are unentangled. By contrast, operators on any row/column 
of MS define an entangled base. A square/grid like the MS was used by Mermin [2] — and frequently 
referred to as a Mermin's square since then — to provide a simple proof of the Kochen-Specker theorem 
in four dimensions. The proof goes as follows. One observes that the square is polarized in the sense that 
the product of three operators on any column equals +I4 (the 4x4 identity matrix), while the product 
of three observables on any row equals — I4. By multiplying all columns and rows one gets — 14. This is, 
however, not the case for the eigenvalues of the observables; they all equal ±1 and their corresponding 
products always yield +1 because each of them appears in the product twice; once as the eigenvalue in 
a column and once as the eigenvalue in a row. The algebraic structure of mutually commuting operators 
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Figure 3: The Mermin square MS viewed as a "sub-Pappus" configuration; the Pappus configuration (93) 
is obtained by adding the three extra lines (dotted). 

thus contradicts that of their eigenvalues, which furnishes a proof of the Kochen-Specker theorem. The 




(PG) (I) 

Figure 4: The partitioning of V[2, 2] into a maximum independent set (/) and the Petersen graph (PG), 
aka its minimum vertex cover. The two vertices on an edge of PG correspond/map to a vertex in / (as 
illustrated by the labels on the edges of a selected closed path) . 

MS set is also recognized as a (92, 63) configuration for any point is incident with two lines and any line 
is incident with three points and does not change its shape if we reverse our notation, i. e., join by an edge 
two mutually non-commuting observables; in graph theoretical terms this means that the MS equals its 
complement. It is also interesting to see that this configuration sits inside the Pappus (93) configuration 
(all vertices and lines in Fig. 3) by removing from the latter the three non-concurrent lines (the dotted 
ones). Last but not least, it needs to be mentioned that the MS configuration represents also the structure 
of the projective line over the product ring Z2 X Z2 if we identify the points sets of the two and regard 
edges as joins of mutually distant points [141 115) : it was precisely this fact that motivated our in-depth 
study of projective ring lines [131 117] and finally led to the discovery of the relevant geometries behind 
two- and multiple-qubit systems [3 El [7] . 

3.3 The Petersen graph PG and the maximum independent set I 

The third fundamental partitioning of V[2,2] comprises a maximum independent set / and the Petersen 
graph PG [5J. This can be done in six different ways and one of them features / = (1,2,6,9,12) and 
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Figure 5: The complement of the Petersen graph viewed as the Desargues configuration; every line com- 
prises three pairwise non-commuting operators Ofe, o;, o m , k ^ I ^ m, i.e., the operators obeying the rule 
Ok-oi = ±io m . 

PG = (3, a, 4, 5, b, 7, 8, c, 10, 11). As in the case of their cousins CB and BP, the Petersen graph PG 
admits a map of its edges into the vertices of the independent set I. Its complement, PG, can be viewed 
as a Desargues configuration (IO3) (see Fig. 5) whose points are the vertices of PG and lines are triples 
of raon-commuting observables Ok, 01, o m , k ^= I ^= m, o^.oi — ±io m . The Desargues configuration is, like 
those of Fano and Pappus, self-dual. 

3.4 Finite projective algebraic geometry underlying V[2,2] 
3.4.1 V[2, 2] as the generalized quadrangle of order two — W(2) 

At this point we have dissected 7^2,2] to such an extent that we are ready to show the unique finite 
projective geometry hidden behind — namely the generalized quadrangle of order two, W(2) [6|. As 
already mentioned in Sec. 2.2, W(2) is the simplest thick generalized quadrangle endowed with fifteen 
points and the same number of lines, where every line features three points and, dually, every point is 
incident with three lines, and where every point is joined by a line (or, simply, collinear) with other 
six points [TSldO]. These properties can easily be grasped from the drawing of this object, dubbed for 
obvious reasons the doily, depicted in Fig. 6; here, all the points are drawn as small circles, while lines 
are represented either by line segments (ten of them), or as segments of circles (the remaining five of 
them). To recognize in this picture V[2, 2] one just needs to identify the fifteen points of W{2) with our 
fifteen generalized Pauli operators as explicitly illustrated, with the understanding that collinear means 
commuting (and, so, raon-collinear reads nora-commuting); the fifteen lines of W(2) thus stand for nothing 
but fifteen maximum subsets of three mutually commuting operators each. 

That W{2) is indeed the right projective setting for V[2, 2] stems also from the fact that it gives a 
nice geometric justification for all the three basic partitionings/factorizations of V[2,2]. To see this, we 
just employ the fact that W(2) features three distinct kinds of geometric hyperplanes [18]: 1) a perp-set 
(H c i(X)), i.e., a set of points collinear with a given point X, the point itself inclusive (there are 15 such 
hyperplanes); 2) a grid (H gr ) of nine points on six lines, aka a slim generalized quadrangle of order (2, 1) 
(there are 10 such hyperplanes); and 3) an ovoid (H ov ), i.e., a set of (five) points that has exactly one point 
in common with every line (there are six such hyperplanes). One then immediately sees that a perp-set 
is identical with a Fano pencil, a grid answers to a Mermin square and, finally, an ovoid corresponds to a 
maximum independent set. Because of self-duality of W(2), each of the above introduced hyperplanes has 
its dual, line-set counterpart. The most interesting of them is the dual of an ovoid, usually called a spread, 
i. e., a set of (five) pairwise disjoint lines that partition the point set; each of six different spreads of W(2) 
represents such a pentad of mutually disjoint maximally commuting subsets of operators whose associated 
bases are mutually unbiased [3j|4j. It is also important to mention a dual grid, i.e., a slim generalized 
quadrangle of order (1,2), having a property that the three operators on any of its nine lines share a base 
of unentagled states. It is straightforward to verify that these lines are defined by the edges of a BP; 
each of the remaining six lines (fully located in the corresponding/complementary MS) carries a base of 
entangled states (see Fig. 6). 
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Figure 6: W(2) as the unique underlying geometry of two-qubit systems. The Pauli operators correspond 
to the points and maximally commuting subsets of them to the lines of the quadrangle. Three operators 
on each line have a common base; six out of fifteen such bases are entangled (the corresponding lines being 
indicated by boldfacing). 



We shall finish this section with the following observation. A triad of a generalized quadrangle is an 
unordered triple of pairwise non-collinear points, with the common elements of the perp-sets of all the 
three points called its centers [18]. W{2) possesses two different kinds of triads: 1) those featuring three 
centers (e. g., the triple {&, 5, 11}), as well as 2) those which are imicentric (e. g., the triple {1, 6, 12}). 

3.4.2 V[2,2] and the projective line over the full two-by-two matrix ring over Z 2 

W(2) is found as a sw&geometry of many interesting projective configurations and spaces [TBI 120] ■ We 
will now briefly examine a couple of such embeddings of 1^(2) in order to reveal further intricacies of its 
structure and, so, to get further insights into the structure of the two-qubit Pauli graph. 

We shall first consider an embedding of 14^(2) in the projective line defined over the ring Z\ x2 of full 
2x2 matrices with Z 2 -valued coefficients, 

Z 2x2 ^{( « J ) | a,M,5EZ 2 ^, (1) 

because it was this projective ring geometrical setting where the relevance of the structure VT(2) for two- 
qubits was discovered [B]. To facilitate our reasonings, we label the matrices of Z 2 x2 in the following 
way 



1 j ' — \ 1 j ' v 1 1 J v 1 1 



1 j ' \ 1 j 7 \o y ' V 

; ; v 10'^ r M,n^n o Vi2'.^ 1 



1 j ' \i j * \ 1 1 j ' \i 1 

;)<"' s (o ;)' i5 ' s (i o)'°' s (o 0). < 2 > 

and see that {1', 2', 9', 11', 12', 13'} are units (i. e., invertible matrices) and {0', 3', 4', 5', 6', 7', 8', 10', 14', 15'} 
are zero-divisors (i.e., matrices with vanishing determinants), with 0' and 1' being, respectively, the ad- 
ditive and multiplicative identities of the ring. Employing the definition of a projective ring line given in 
Sec. 2.2, it is a routine, though a bit cumbersome, taslQ to find out that the line over Z 2x2 is endowed 



1 See, for example, |13l 1141 for more details about this methodology and a number of illustrative examples of a projective 
ring line. 
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Figure 7: A BP + MS factorization of V[2,2}) in terms of the points of the subconfiguration of the 
projective line over the full matrix ring Z 2x2 ; the points of the BP have both coordinates units, whilst 
those of the MS feature in both entries zero-divisors. The "polarization" of the Mermin square is in this 
particular ring geometrical setting expressed by the fact that each column/row is characterized by the 
fixed value of the the first/second coordinate. Compare with Fig. 2. 



with 35 points whose coordinates, up to left-proportionality by a unit, read as follows 
(l',l'), (l',2'), (l',9'), (l',ll'), (l',12'),(l',13'), 

(l',0'), (l',3'), (l',4'), (l',5'), (l',6'), (l',n (l',8'), (IMC), (l',14'), (1', 15'), 
(0',l'), (3',1'), (4',1'), (5M'), (6',1'), (7U'), (8M'), (10', 1'), (W, 1'), (15', 1'), 
(3', 4'), (3', 10'), (3', 14'), (5', 4'), (5', 10'), (5', 14'), (6', 4'), (6', 10'), (6', 14'). (3) 

Next, we pick up two mutually distant points of the line. Given the fact that GL(2, R) acts transitively 
on triples of pairwise distant points [16] , the two points can, without any loss of generality, be taken to 
be the points Uq := (1, 0) and Vb := (0, 1). The points of W(2) are then those points of the line which are 
either simultaneously distant or simultaneously neighbor to Uo and Vo- The shared distant points are, in 
this particular representation, (all the) six points whose both entries are units, 

(l',l'), (l',2'), (l',9'), 

(l',ll'), (l',12'), (l',13'), (4) 

whereas the common neighbors comprise (all the) nine points with both coordinates being zero-divisors, 

(3', 4'), (3', 10'), (3', 14'), 
(5', 4'), (5', 10'), (5', 14'), 

(6', 4'), (6', 10'), (6', 14'), (5) 

the two sets thus readily providing a ring geometrical explanation for a BP + MS factorization of the 
algebra of the two-qubit Pauli operators, Fig. 7, after the concept of mutually neighbor is made synonymous 
with that of mutually commuting [6]. To see all the three factorizations within this setting it suffices to 
notice that the ring Z 2x2 contains as subrings all the three distinct kinds of rings of order four and 
characteristic two, viz. the (Galois) field F4, the local ring Zi\x\j (x 2 ), and the direct product ring Z2 x Z2 
[21j , and check that the corresponding lines can be identified with the three kinds of geometric hyperplancs 
of W{2) as shown in Table 5 0. 

The other embedding of W{2) to be briefly dealt with is the one into the projective space, PG(3, 2), as 
illustrated in Fig. 8. This embedding is, in fact, a very close ally of the previous one due to a remarkable 
bijective correspondence between the points of the line over Z 2x2 and the lines of PG(3,2) [52]. W{2) 
and PG(3, 2) are identical as the point sets, whilst the fifteen lines of W(2) are so-called totally isotropic 
lines with respect to a symplectic polarity of PG(3, 2) (Sec. 4.2). 
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Figure 8: An illustration of an embedding of the generalized quadrangle W(2) (and thus of the associated 
Pauli graph V[2, 2]) into the projective space PG(3, 2). The points of PG(3, 2) are the four vertices of the 
tetrahedron, its center, the four centers of its faces and the six centers of its edges; the lines are the six 
edges of the tetrahedron, the twelve medians of its faces, the four circles inscribed in the faces, the three 
segements linking opposite edges of the tetrahedron, the four medians of the terahedron and, finally, six 
circles located inside the tetrahedron [20 . The fifteen points of PG(3, 2) correspond to the fifteen Pauli 
operators/ vertices of V[2, 2]. All the thirty-five lines of the space carry each a triple of operators Ofe, oi, 
o m , k 7^ I 7^ m, obeying the rule o^.o; = no m ; the operators located on the fifteen totally isotropic lines 
belonging to W(2) yield fj, = ±1, whereas those carried by the remaining twenty lines (not all of them 
shown) give /i = ±i. 

4 The Pauli graph of iV-qubits 

Following the same strategy as in the preceding section, we find out that the 4 3 — 1 = 63 tensor products 
Cj <8> crj (g) crfc, [i, j, k = 1, 2, 3, 4, (i, j, k) ^ (1, 1, 1)] form the vertices and their commuting pairs the edges of 
a regular graph of degree 30, P[2,3], with spectrum {— 5 27 , 3 35 , 30}. The corresponding incidence matrix 
can also be cast into a compact tripartite form, Table 6, after the reference points 03 = a x ® I2 <8> I2, 
&3 = G y ® I2 <£> I2 and C3 = a z (E> I2 ® I2 have been omitted. This matrix looks very much the same as its 
two-qubit counterpart (Table 3), save for the fact that now all the submatrices are of rank 15 x 15. As 
in the two-qubit case, the matrix ^3 can simply be viewed as the join of O3 and the unit matrix Ig. The 
same self-similarity pattern interrelating the incidence matrices of (N + 1)- and A^-qubit systems is found 
for any N. 

As for the two-qubit incidence matrix, one of the most natural factorizations of the three-qubit matrix 
consists of the first block O3 and a larger square block M3, of cardinality 45, containing O3 and A3. The 
latter block is self-complementary, as is its two-qubit counterpart, a Mermin square; it represents a regular 
graph of degree 22 and spectrum {— 5 10 , — 3 9 , — 2 2 , l 5 , 3 18 , 22}. The structure of this block is very intricate: 



Table 5: Three kinds of the distinguished subsets of the generalized Pauli operators of two-qubits (V[2, 2])) 
viewed either as the geometric hyperplanes in the generalized quadrangle of order two (W(2)) or as the 
projective lines over the rings of order four and characteristic two residing in the projective line over Z 2 x2 . 



T[2,2] 


set of five mutually 


set of six operators 


nine operators of a 




non-commuting operators 


commuting with a given one 


Mermin's square 


W(2) 


ovoid 


perp-set\ {reference point} 

Z2[X}/(X 2 ) 


grid 


Proj. Lines over 


F± = Z 2 [x]/(x 2 +x + l) 


Z 2 x Z 2 = Z 2 [x]/(x(x + 1)) 
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it can be recovered again by removing from the reduced incidence matrix shown in Table 6 the first triple 
of points and all the reference points (of the type a, b and c, see Table 2) of the parent scale, i. e., an extra 
set of 3 + 3 x 4 = 15 "pseudo-reference" points of the "daughter" scale. 



o 3 


A 3 


A 3 


A 3 


A 3 


o 3 


A 3 


A 3 


A 3 


A 3 


o 3 


A 3 


A 3 


A 3 


A 3 


o 3 



Table 6: The incidence matrix of "P[2,3] after removal of the triple of reference points (compare with 
Table 3). 



After a closer look at A^3, one reveals in it three subsets isomorphic to the Mermin square of two-qubits 
(Fig. 2 and/or Fig. 7), from which we can form doubles (18 points) and triples (27 points) having spectra 
{— 3 4 , — l 8 , 0, 3 4 , 8} and {— 3 12 , 6 , 3 8 , 12}, respectively; the graph of the latter bears number 105 in the list 
of graphs with few eigenvalues given in [23] . One can also form m-tuples of the "generalized" Mermin square 
of size m — 1,2,3,4 using the "entangled" subset £ located in the first block 3 and the extra M S copies 
from A4 3 , to get another interesting blocks SUMS, £U(2x MS) and £U(3x MS) and the associated graphs 
with spectra {-3 4 , -l 9 , 3 4 , 9}, {-3 12 , 5 , 3 8 , 3(2 ± v^)} and {-5 4 , -3 12 , 2 , l 4 , 3 12 , 8 ± y/§l}, respectively. 

4.1 Rank N symplectic polar spaces behind the iV-qubit Pauli graphs 

The geometry underlying higher order qubits [7] can readily be hinted from the observation that our doily 
W(2), embodying the two-qubit operators' algebra, is the lowest rank representative of a big family of 
symplectic polar spaces of order two. 

A symplectic polar space (see, e.g., [THJ[2I1[25] for more details) is a d-dimensional vector space over 
a finite field F q , V(d, q), carrying a non-degenerate bilinear alternating form. Such a polar space, usually 
denoted as Wd-i{q), exists only if d = 2N, with TV being its rank. A subspace of V(d, q) is called totally 
isotropic if the form vanishes identically on it. W2n-i{q) can then be regarded as the space of totally 
isotropic subspaces of PG(2N — 1, q) with respect to a symplectic form, with its maximal totally isotropic 
subspaces, called also generators G, having dimension N — 1. For q = 2, this polar space contains 

|WW-i(2)| = \PG{2N - 1,2)| = 2 2N - 1 = 4 N - 1 (6) 

points and (2 + 1)(2 2 + 1) . . . (2 N + 1) generators. A spread S of W2n-i(q) is a set of generators partitioning 
its points. The cardinalities of a spread and a generator of W^iV-i(2) read 

|S|=2 W + 1 (7) 

and 

\G\ =2 N -1, (8) 

respectively. Finally, it needs to be mentioned that two distinct points of W^at-i (q) are called perpendicular 
if they are joined by a line; for q = 2, there exist 

#A = 2 2N - 1 (9) 

points that are not perpendicular to a given point. 

Now, in light of Eq. (6), we can identify the Pauli operators of iV-qubits with the points of H^jv-i(2). If, 
further, we identify the operational concept "commuting" with the geometrical one "perpendicular," from 
Eqs. (7) and (8) we readily see that the points lying on generators of W2jv-i(2) correspond to maximally 
commuting subsets (MCSs) of operators and a spread of W2jv-i(2) is nothing but a partition of the whole 
set of operators into MCSs. Finally, Eq. (9) tells us that there are 2 2N ~ 1 operators that do not commute 
with a given operator^ 

Recognizing W / 2Ar-i(2) as the geometry behind iV-qubits, we will now turn our attention on the prop- 
erties of the associated Pauli graphs, V[2, N]. 



2 Shortly after Ref. [7] was posted on the arXiv-e, physicist D. Gross (Imperial College, London) sent us an outline of the 
proof of this property and a couple of weeks later, Koen Thas (Ghent University), a young mathematician, also informed us 
about finding a proof of the same statement. 
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4.2 Strong regularity of the iV-qubit Pauli graph 



As already introduced in Sec. 2.1, a strongly regular graph, srg(i;, D, A, /i), is a regular graph having v 
vertices and degree D such that any two adjacent vertices are both adjacent to a constant number A of 
vertices, and any two distinct non-adjacent vertices are also both adjacent to a constant number fi of 
vertices. It is known that the adjacency matrix A of any such graph satisfies the following equations [26j 



AJ = DJ, A 2 + {n - X)A + D)I = fiJ, 



(10) 



where J is the all-one matrix. Hence, A has D as an eigenvalue with multiplicity one and its other 
eigenvalues are r (> 0) and I (< 0), related to each other as follows: r + I = A — fx and rl = fx — D. 
Strongly regular graphs exhibit many interesting properties [26 . In particular, the two eigenvalues r and 
I are, except for (so-called) conference graphs, both integers, with the following multiplicities 



-D(l + l)(D-l) _ D{r + l){D-r) 

1 {D + rl)(r-l) g {D + rl)(r-l) : 



(11) 



respectively. The A-qubit Pauli graph is strongly regular, and its properties can be inferred from the 
relation between symplectic polar spaces and partial geometries. 

A partial geometry is a more general object than a finite generalized quadrangle. It is finite near-linear 
space {P, L} such that for any point P not on a line L, (i) the number of points of L joined to P by a 
line equals a, (ii) each line has (s + 1) points, (iii) each point is on (t + 1) lines; this partial geometry 
is usually denoted as pg(s,i, a) [11]. The graph of pg(s,i, a) is endowed with v — (s + 1) ^'^"^ vertices, 
C = (t + 1) i st+a ) lines and is strongly regular of the type 



srg (s + 1) 



(st + a) 



a 



, s(t + 1), s - 1 + t(a - l),a(t + 1) 



(12) 



The other way round, if a strongly regular graph exhibits the spectrum of a partial geometry, such a 
graph is called a pseudo-geometric graph. Graphs associated with symplectic polar spaces W2n-i(q) are 
pseudo-geometric [35], being 



pg q 



q 



N-l 



1 



g-i 



• q 



N-l 



q 



N-l 



1 



q-l 



-graphs. 



(13) 



Combining these facts with the findings of the preceding section, we conclude that that 7V-qubit Pauli 
graph is of the type given by Eq. (13) for q — 2; its basics invariants for a few small values of N are listed 
in Table 7. 



N 


V 


C 


D 


r 


I 


A 


M 


s 


t 


a 


2 


15 


15 


6 


1 


-3 


1 


3 


2 


2 


1 


3 


63 


45 


30 


3 


-5 


13 


15 


6 


4 


3 


4 


255 


153 


126 


7 


-9 


61 


63 


14 


8 


7 



Table 7: Invariants of the Pauli graph "P[2, N], N = 2, 3 and 4, as inferred from the properties of the 
symplectic polar spaces of order two and rank N. In general, i 



iJV 



1, D 



1-2 



2N-1 



2-1 



■)N-1 



t 

can also be found from s, t and a themselves. 



[i = a(t + 1) = rl + D and A = s — 1 + t(a — 1)) = fi + r + 1. The integers v and e 



5 The Pauli graph of two-qutrits 

A complete orthonormal set of operators of a single-qutrit Hilbert space is [5] 

G! = {I 3 ,Z,X,Y,V,Z 2 ,X 2 ,Y 2 ,V 2 }, 1= 1,2, 3,..., 9, (14) 

/ 1 \ / 1 \ 

where 7 3 is the 3 x 3 unit matrix, Z = w ,A=10 ,F = XZ, V = XZ 2 and 

\ uj 2 J \ 1 ) 

oj = cxp(2i7r/3). Labelling the two-qutrit Pauli operators as follows 1 = I3 ® a\, 2 = I3 ® (T2, •"> 
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Figure 9: A partitioning of Wg into a grid (top left), an 8-coclique (top right) and a four-dimensional 
hypercube (bottom). 
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Figure 10: A partitioning of Wg into a tripartite graph comprising a 10-coclique, two 9-cocliques and a 
set of four triangles; the lines corresponding to the vertices of a selected triangle intersect at the same 
observables of Vg and the union of the latter form a line of Vg. 

8 = I3 <g> <7 8 , a = <7i <g) I 3 , 9 = ai ® . . , b = <t 2 ® I3, 17 = a 2 ® cr lv . . , c = <t 3 ® I3,. . ., h = <r 8 (8) 7 2) - ■ 
72 = er 8 (g) erg , one obtains the incidence matrix of the two-qutrit Pauli graph Vg. 

Computing the spectrum {— 7 15 , — l 40 , 5 24 , 25} one observes that the graph is regular, of degree 25, 
but not strongly regular. The structure of observables in Vg is much more involved than in the case of 
two-qubits although it is still possible to recognize identifiable regular subgraphs. In order to get necessary 
hints for the geometry behind this system, it necessitates to pass to its dual graph, Wg, i.e., the graph 
whose vertices are maximally commuting subsets (MCSs) of Vg. To this end, let us first give a complete 
list of the latter: 

Ll = {1, 5, a, 9, 13, e, 41, 45}, L 2 = {2, 6, a, 10, 14, e, 42, 46}, L 3 = {3, 7, a, 11, 15, e, 43, 47}, 
L 4 = {4, 8, a, 12, 16, e, 44, 48}, Mi = {1, 5, b, 17, 21, /, 49, 53}, M 2 = {2, 6, b, 18, 22, /, 50, 54}, 
M 3 = {3, 7, b, 19, 23, /, 51, 55}, M 4 = {4, 8, b, 20, 24, /, 52, 56}, JVi = {1, 5, c, 25, 29, g, 57, 61}, 
N 2 = {2, 6, c, 26, 30, g, 58, 62}, N 3 = {3, 7, c, 27, 31, g, 59, 63}, N 4 = {4, 8, c, 28, 32, g, 60, 64}, 
Pi = {1, 5, d, 33, 37, h, 65, 69}, P 2 = {2, 6, d, 34, 38, h, 66, 70}, P 3 = {3, 7, d, 35, 39, h, 67, 71}, 

P 4 = {4, 8, d, 36, 40, h, 68, 72}, 
Xi = {9, 22, 32, 39, 45, 50, 60, 67}, X 2 = {10, 17, 27, 40, 46, 53, 63, 68}, X 3 = {11, 20, 30, 33, 47, 56, 58, 69}, 
X 4 = {12, 23, 25, 34, 48, 51, 61, 70}, X 5 = {13, 18, 28, 35, 41, 54, 64, 71}, X 6 = {14, 21, 31, 36, 42, 49, 59, 72}, 

X 7 = {15, 24, 26, 37, 43, 52, 62, 65}, X 8 = {16, 19, 29, 38, 44, 55, 57, 66}, 
Yi = {9, 23, 30, 40, 45, 51, 58, 68}, Y 2 = {10, 19, 32, 33, 46, 55, 60, 69}, Y 3 = {11, 22, 25, 36, 47, 50, 61, 72}, 
Ya = {12, 17, 26, 39, 48, 53, 62, 67}, Y 5 = {13, 20, 27, 34, 41, 56, 63, 70}, Y 6 = {14, 23, 28, 37, 42, 51, 64, 65}, 



14 




Y T = {15, 18, 29, 40, 43, 54, 57, 68}, Y a = {16, 21, 30, 35, 44, 49, 58, 71}, 
Zi = {9, 24, 31, 38, 45, 52, 59, 66}, Z 2 = {10, 24, 25, 35, 46, 52, 61, 71}, Z 3 = {11, 17, 28, 38, 47, 53, 64, 66}, 
Z 4 = {12, 18, 31, 33, 48, 54, 59, 69}, Z 5 = {13, 19, 26, 36, 41, 55, 62, 72}, Z 6 = {14, 20, 29, 39, 42, 56, 57, 67}, 
Z 7 = {15, 21, 32, 34, 43, 49, 60, 70}, Z s = {16, 22, 27, 37, 44, 50, 63, 65}. 

From there we find that Wg consists of 40 vertices and has spectrum {— 4 15 , 2 24 , 12}, which are the char- 
acteristics identical with those of the generalized quadrangle of order three formed by the totally singular 
points and lines of a parabolic quadric 0(4,3) in PG(4, 3)[TH]- The quadrangle Q(4, 3), like its two-qubit 
counterpart, exhibits all the three kinds of geometric hyperplanes, viz. a slim generalized quadrangle of 
order (3,1) (a grid), an ovoid, and a perp-set, and these three kinds of subsets can all indeed be found to 
sit inside Wg. One of the grids is formed by the sixteen lines Li, Mi, Ni and Pi (i — 1, 2, 3 and 4) as 
illustrated in Fig. 9; the remaining 24 vertices comprise an 8-coclique (Xi, which correspond to mutually 
unbiased bases), and a four-dimensional hypercube {Yi and Zi). Next, one can partition Wg into a maxi- 
mum independent set and the minimum vertex cover using a standard graph software. The cardinality of 
any maximum independent set is 10 (= 3 2 + 1), which means that any such set is an ovoid of 0(4, 3) [15]. 
It is easy to verify that, for example, the set {Li, M 2 , N 3 , P4, X 3 , X 8 , 14, Y 6 , Z 2 , Z7} is an ovoid; given any 
maximum independent set, Q{4, 3)/Wg can be partitioned as shown in Fig. 10. The remaining type of a 
hyperplane of Wg is a perp-set, i.e. the set of 12 vertices adjacent to a given ("reference") vertex (Fig. 11); 
the set of the remaining 27 vertices can be shown to consist of three ovoids which share (altogether and 
pairwise) just a single vertex — the reference vertex itself. This configuration bears number 99 in a list of 
graphs with few eigenvalues given in Ref. [23] and can schematically be illustrated in form of a "triangle" , 
with a triangular pattern at its nodes and a 1 x 2 grid put on its edges; the union of 1 x 2 grid and a 
triangle cither forms a Mermin-square-type graph M , as already encountered in the two-qubit case, or a 
quartic graph of another type, denoted as K (see Fig. 11). 

The foregoing observations and facts provide a reliable basis for us to surmise that the geometry 
behind Wg is identical with that of Q(4,3). If this is so, then the symplectic generalized quadrangle of 
order three, W(3), which is the dual of Q(4, 3)[TH], must underlie the geometry of the Pauli graph Pg. 
However, the vertex-cardinality of W(3) is 40 (the same as that of Q(4, 3)), whilst Pg features as many 
as 80 points/vertices. Hence, if the geometries of W(3) and Pg are isomorphic, then there must exits a 
natural pairing between the Pauli operators such that there exists a bijection between pairs of operators of 
Pg and points of W(3). This issue requires, obviously, a much more elaborate analysis, to be the subject 
of a separate paper [34], 



6 Conclusion 

The paper introduces an important concept of the Pauli graph for the generalized Pauli operators of 
finite-dimensional quantum systems and illustrates and discussed this concept in an exhaustive detail for 
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iV-qubit systems, N > 2. In doing so, the geometries underlying these systems, viz. the symplectic polar 
spaces of rank TV and order two, are invoked to reveal all the intricacies of the algebra of the operators and 
its basic factorizations. Although there exits a variety of other interesting geometry-oriented approaches 
to model finite dimensional quantum systems (see, for example, 27J [33 ), ours seems to be novel in that it 
goes beyond classical projective geometry and Galois fields and is, in principle, applicable to any quantum 
system of finite dimension. 
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